Droplets and the Configurational Entropy Crisis for Random First Order Transitions 
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We consider the effect of droplet excitations in the random first order transition theory of glasses 
on the configurational entropy. The contribution of these excitations is estimated both at and 
above the ideal glass transition temperature. The temperature range where such excitations could 
conceivably modify or "round-out" an underlying glass transition temperature is estimated, and 
found to depend strongly on the surface tension between locally metastable phases in the supercooled 
liquid. For real structural glasses this temperature range is found to be very narrow, consistent with 
the quantitative success of the theory. For certain finite-range spin-glass models, however, the surface 
tension is estimated to be significantly lower leading to much stronger entropy renormalizations, thus 
providing an explanation for the lack of a strict thermodynamic glass transition in simulations of 
these models. 



A number of models of disordered glassy systems such 
as random heteropolymers pi , p-spin glasses M and Potts 
spin glasses f| in the mean field limit undergo a ther- 
modynamic transition at which an appropriately defined 
configurational entropy vanishes. While supercooled liq- 
uids in the laboratory never reach such a transition be- 
fore falling out of equilibrium, rather natural extrapo- 
lations of the experimentally measured thermodynamics 
also suggest that if a fluid could remain equilibrated upon 
slow cooling, an entropy crisis would occur at an ideal 
glass transition temperature Tk B . This temperature is 
remarkably near the temperature where, again a natural 
extrapolation of the experimental kinetic data suggests 
that relaxations become infinitely slow. 

Over many decades, contemplating such an ideal glass 
transition has entered into discussions of glassy dynam- 
ics. The "random first order transition" (RFOT) exhib- 
ited by the mean-field models has become a central ele- 
ment of the modern theory of glasses ||, 0, [|, 0, fo|. 
Recently this theory has been shown to predict quan- 
titatively from a microscopic starting point the dynam- 
ics of supercooled liquids in the laboratory |pd| , |l2|| , in- 
cluding the relaxation time distribution and the length 
scale of dynamical heterogeneity Jl3). This theory also 
leads to a natural way to describe the anomalous ther- 
modynamic and transport properties of low temperature 
glasses, when the theory is quantized |T^]. Despite its 
venerable history |lq] the idea of an underlying ideal 
glass transition related to an entropy crisis has also been 
greeted by much controversy since it requires extrapola- 
tion from extant data. Yet only one a priori objection to 
an entropy crisis has been crisply formulated: configura- 
tional entropy cannot vanish above absolute zero since a 
system with finite range forces can always support point 
defects at any finite temperature |T^| . The most naive 
incarnation of this objection, the existence of vacancies 
or interstitials, can easily be dismissed by practical per- 
sons. At the densities of glasses the concentration of 
such defects is miniscule. Reasonable estimates suggest 
they occur in f part in 10 4 so they hardly contribute to 
the configurational entropy, which is typically Iks per 



particle at the laboratory liquid-glass transition temper- 
ature, T g . More troubling is that the random first order 
transition theory itself requires the existence of excita- 
tions ("entropic droplets") which are localized and have 
finite energy cost, ft is conceivable these excitations can 
smear out or modify the critical exponents of the tran- 
sition predicted by mean-field theory. In this paper, we 
estimate what the droplet contribution to the entropy 
both below and above an ideal glass transition would be. 
This calculation parallels Fisher's droplet calculations on 
the essential singularities in the thermodynamics encoun- 
tered near an ordinary first order phase transition [JL7|. 
The calculation shows how close to an ideal glass transi- 
tion we must go to see these renormalization effects and 
quantifies where the notion of configurational entropy be- 
comes ambiguous. In essence this calculation provides 
something like a "Ginzburg" criterion jL8| describing the 
breakdown of mean-field theory for a random first order 
transition. The size of the strongly renormalized regime 
depends strongly on the effective surface tension between 
mean-field replica symmetry broken phases. The effective 
surface tension parameter is actually quite large for the 
liquid-glass transition because it depends on the Linde- 
mann vibrational amplitude which is so small. Effectively 
there is, then, a small parameter in the theory which al- 
lows the mean-field estimates to prevail until very close to 
the transition (much as does the BCS theory of conven- 
tional superconductivity). The entropy crisis for liquids 
will be noticeably avoided only when the residual entropy 
is of the order 10 _4 A:b and (T — Tk)/Tk is approximately 
10 fee/Ac,, where Ac v is the mean-field heat capacity 
discontinuity at the transition. In contrast, the effective 
surface tension for a 10-state nearest neighbour lattice 
Potts glass is quite small. This small value is consistent 
with the avoidance of an entropy crisis in a recent simu- 
lation of this system [|l9). Since the surface tension pa- 
rameter also tunes the degree of slowing as the transition 
is approached understanding these effects may allow the 
design of simulations that can more effectively confront 
the theory. 

We begin by very briefly summarizing the essential 
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points of the RFOT theory of glasses, and then construct 
a simple partition function for the droplet excitations. 
After discussion of droplet effects at Ik, we evaluate the 
temperature range over which such excitations will be sig- 
nificant for a structural glass as T approaches Xk from 
above. Finally, we apply the results to several glassy 
systems with different surface tensions. 

Below a temperature Ta, which corresponds to the 
dynamical glass transition predicted by mode-coupling 
theory |po| , reconfiguration of a real supercooled liquid 
occurs by activated processes. According to the RFOT 
theory of glasses for Tk < T < Ta the driving force for 
reconfiguration of a locally frozen region is the configura- 
tional entropy, while there is an energetic cost due to the 
surface tension between two different locally metastable 
solutions. The free energy profile for reconfiguration of a 
region of radius r is thus 

F(r) = 4:irr 2 a{r) - ^nr 3 pTS c (1) 

where S c is the configurational entropy density per 
molecule, and p is the number density. For the ex- 
perimental comparison the measured S c was used, but 
the theory more formally would call for its mean-field 
value. Generalization of a renormalization group (RG) 
argument due to Villain ]2l| suggests that the surface 
tension is reduced from its bare value <jq (obtained be- 
tween two typical metastable states) due to wetting of 
the interface by additional metastable states, accord- 
ing to a(r) — CTo( r o/ r ) 5 ; with ro — a typical in- 
termolecular distance ||. It is helpful to express the 
two key parameters (surface tension and configurational 
entropy) in terms of dimcnsionlcss parameters that are 
approximately constant as T — > Tk- Specifically, the 
configurational entropy is assumed to vanish linearly, 
S c = T ~J K -j^-, where Ac v is the heat capacity jump 
per particle at the transition, and we define a reduced 
surface tension A = 47rrgcro(TK)/fcBTk. In the region 
T K < T < Ta, we assume that 47rrgcro(r)/fc s T is slowly 
changing, and may be approximated by A. The recon- 
figuration barrier may then be written in Vogel-Fulcher 
form as /3AF* = DT K /(T - T K ), with the fragility 
D = Ae j kB ■ The supercooled liquid is viewed as 

a mosaic of metastable regions of typical size R/ro = 
— - 

(ffsf;) 3 (r^zk) 3 ' as follows irom F ( R ) = For 
structural glasses Ac„/fcB varies a great deal from sub- 
stance to substance over a typical range of 0.2 to 5. How- 
ever, a density functional calculation [|llj gives A « 23 for 
structural glasses approximately independently of molec- 
ular details since it depends only logarithmically on the 
Lindemann ratio of mean-square vibrational amplitude 
to near-neighbour spacing which characterizes the vibra- 
tions at T g . 

Within each domain of the mosaic there should ex- 
ist further droplet excitations with free energy profiles 
also given by eq. (|lj), as shown schematically in fig. 
fil. Indeed, since these droplets may also contain even 




FIG. 1: A schematic of the structure of a supercooled liquid 
near Tk, showing droplet excitations within the mosaic cells. 



smaller droplets, this suggests that like the surface ten- 
sion the configurational entropy should be determined 
by RG equations. In this paper, we estimate the re- 
gion where such a renormalization should be significant 
and leave the full RG calculation for future work, al- 
though we note that in the case of the random field 
Ising model, the presence of droplets within droplets 
is not sufficient to raise the lower critical dimension 
from d — 2 to d — 3 |23|. The partition func- 
tion for droplet formation within a single mosaic cell 
is Z = J2{n s } exp(S , com6 ({n i }) - (3 J2i "-i-F(^)), where 
the free energy cost of forming a droplet of radius Ti 
about a given position is given by eq. ([!]), and S com b 
is an entropy associated with the different positions the 
droplets may take up within the domain. Assuming that 
there is approximately one distinct droplet position for 
each molecule that a droplet may be centered on, gives 
—\{pVeff) n * states of rii identical independent droplets. 
The effective volume within a domain in which a droplet 
may move, V e // = ^ttR^j, is taken to be smaller than 
the domain volume due to the presence of the domain 
walls of finite thickness. While strictly speaking the RG 
wetting argument suggests R e f f oc R |8|] , with essentially 
no effect on the results we take R e ff = R — »"o- Ex- 
cluded volume effects have been entirely ignored, even 
between identical droplets, as has the reduction in the 
volume available to position larger droplets due to the 
restriction that the entire droplet must lie within the do- 
main walls (this latter effect may be straightforwardly 
included, but is not significant for parameters where ig- 
noring the excluded volume is a good approximation). 
Under these approximations the entropy S CO mb reduces 
to S comb = J2i { n i m L°K//] - Inn,!}, and the partition 
function may easily be evaluated giving a free energy den- 
sity of independent droplets Fd r0 p = — {ksT/pV) In Z = 
—kT(V e ff/V)J2i e ~ l3Fi - Droplets may be distinguished 
by the number of molecules they contain, so what is 
meant by the sum over droplet types i in eq. (^) is 
i = 1, 2, . . . , i max , where i max is the maximum size of 
a droplet that can fit inside the domain i max < pV e ff] 
alternatively — and no more approximately — we may 
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make the replacement Y^l 



•ff „2 



r • dr, with 



T s = (Jj)^ r o ~ 0.62ro. The corresponding droplet en- 
tropy density (which follows from differentiation with S c 
and <t fixed) is 



Sdrop/ks = 47TP 



eff 



V 



r 2 {l + A{r/r )i)e- pF{r) dr. (2) 



At Tk the droplet entropy density Sdrop^K) /&b = 



(87T/3) 



3f|/A + 3/A 2 exp(-Af s 2 ), where 



r s /rQ. For a structural glass, with the values given above, 
Sdro P {T}^)/kB ~ 3 x 10~ 5 . For structural glasses, the con- 
tribution to the droplet entropy comes overwhelmingly 
from the smallest droplets and is thus a defect-like en- 
tropy. Unfortunately, on the short length scales of such 
droplets (perhaps involving only a single molecule) the 
calculation of the energetic cost in terms of a surface 
tension becomes less accurate and the details of the par- 
ticular molecules may become relatively more important. 
The effect of the uncertainty in the energetic cost of form- 
ing small droplets on the droplet entropy may be allowed 
for by permitting f s to vary in the range 0.4 to 1.0, which 
corresponds to the energy of the smallest droplets being 
a factor of two smaller or larger than estimated using 
f s = 0.62. Even allowing for this uncertainty the resul- 
tant droplet entropy at Tk remains tiny, in the range 
1 x 10~ 9 to 3 x 10~ 3 k R . 



to 6 x 10~ 6 (for fragile liquids). Such cross-over tem- 
peratures lie well below those accessible to equilibrium 
measurements, since the a relaxation times at T x will be 
of the order To exp(10 6 ), effectively infinite. We show in 
fig. ^| the cross-over temperature as a function of Ac v 
over a range of Ac v appropriate to a structural glass. 
The widening of the renormalized region with decreasing 
Ac v (increasing fragility) is clear. The conclusion that 
T x is inaccessible to direct experimental study for struc- 
tural glasses appears robust to the uncertainties in the 
energy cost calculation: even taking f s — 0.38 (corre- 
sponding to an energy cost of half the natural estimate) 
leads to (Ac v /k B ) x (T x - T K )/T K = 2 x 10~ 3 while 
{Acy/ks) x (T g - T K )/T K » 1. 
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FIG. 2: The configurational entropy and the droplet cor- 
rection to the entropy for a structural glass as a function 
of temperature in the temperature range (T — Tk)/Tk = 
10" 6 fc s /Ac„ to k B /Ac v . 

We now consider T > Tk and estimate the crossover 
temperature, T x , where the droplet entropy and configu- 
rational entropy become comparable. For Tk < T < T x 
droplet effects may modify the character of the random 
first order transition. As SdropiT^) is finite and small, 
it is clear that T x > Tk exists and since, as fig. || shows, 
SdropiT) is not rapidly varying as T approaches Tk, T x 
may be estimated from from SdropiTx) — S C (T). For 
large values of the reduced surface tension and f s = 
0.62 this leads to log 10 [(Ac v /k B ) x (T x - T K )/T K ] = 
—0.21A + 0.3. For structural glasses this suggests (T x — 
Tk)/Tk lies in the range 1.5 x 10~ 4 (for strong liquids) 



FIG. 3: log 10 [(T x - T k )/T k ], where T x is the crossover tem- 
perature as a function of the heat capacity jump, for A = 23 
(bottom; appropriate for a structural glass), A = 12 (middle; 
as may be appropriate for a protein-like random heteropoly- 
mer), and A = 6 (top; possibly appropriate for some spin 
models as discussed in the text). 

While structural glasses all have A w 23 other glassy 
systems are expected to have different bare surface ten- 
sions. The width of the renormalized region depends 
strongly on A. Figure ^ gives the full dependence of 
the width of the renormalized region on the bare sur- 
face tension parameter A. Also shown is the position of 
the laboratory glass transition temperature T g , defined 
in terms of the experimental and microscopic times via 
Tg/T = cxp{DT K /(T g - Tk)) with both r s /r = 10 17 
(applicable to experiments on structural glasses with re- 
laxations in the hours range) and t 9 /tq = 10 s (typical 
of a computer simulation of a spin glass). For large val- 
ues of the surface tension (A > 15) T x is found to lie 
significantly below T g and the additional droplet entropy 
is small at any temperature where equilibration can be 
achieved in practice. For surface tensions in the range 
5 < A < 10, on the other hand, the values of T x and T g 
will be more comparable and the effects of droplet excita- 
tions should become apparent. Depending on the details 
of the model and the timescale accessible to the exper- 
imenter as A is decreased from 10 to 5, at some point 
T x and T g will cross and the effect of the droplet renor- 
malization may be large. We note that by A = 6 the 
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independent droplet approximation breaks down since 
~ 30% of the volume of domains is predicted to be filled 
with droplet excitations at T x . Nevertheless, this does 
not change the conclusion that droplets are significant at 
this surface tension. 

To use the droplet entropy renormalization criterion 
we need to know A for systems with finite-ranged in- 
teractions that undergo RFOTs in the mean-field limit. 
Unfortunately, for protein, Potts and p-spin models accu- 
rate calculations of A are somewhat involved. Therefore, 
here we content ourselves to estimate the value of A for 
various models using rough calculations, that should be 
sufficiently accurate to determine whether or not droplet 
excitations are likely to play an important role in the en- 
tropy. We first make an estimate for polymeric, and then 
some spin-glass models. 




FIG. 4: Logarithm of the width of the renormalized region 
as a function of the dimensionless surface tension parameter. 
The three solid curves correspond (from top to bottom) to 
f s = 0.38, 0.62 and 0.98, with T x calculated from S C (T X ) = 
Sdro P (Tx) and eq. (Q). The three dotted lines are the asymp- 
totes In [(Ac v /k B ) x (T x - T K )/T K ] = -r j A + In [8tt^/3] . 
The two dashed lines give the laboratory glass transition tem- 
perature as log 10 [(Ac„/fc_e) x (T g — Tk)/2k] for waiting times 
appropriate to experiments on structural glasses (lower) and 
computer simulations on spin-glasses (upper). 

The surface tension for a random heteropolymer may 
be estimated using the argument of Takada and Wolynes 
p4| , or alternatively by approximating the surface ten- 
sion between two neighboring misfolded regions of a ran- 
dom heteropolymer with the surface tension of a sharp 
interface between folded and unfolded regions of a pro- 
tein. The latter quantity is a key quantity in the cap- 
illarity model of protein folding p5| , p(| , which in its 
simplest form assumes folding occurs by nucleation with 
a sharp interface between completely folded and com- 
pletely unfolded regions. At the folding temperature, 
X/, the capillarity model predicts a free energy barrier 

of AF/k B T f = ± ' ANi . A knowledge of the vari- 
ation of folding rates at Tf with protein size will thus 
give an estimate of the reduced surface tension. Such 
a study has been performed by Koga and Takada p7j , 
who find for simulations of model proteins a folding rate 



that varies as log 10 kf w — O.ScNs , with c a roughly N- 
independent constant that depends on protein topology 
called relative contact order, c measures how frequently 
the chain returns to itself. The same relation holds for 
real proteins but folding times vary over more than four 
more orders of magnitude and the coefficient is large, 
log 10 kf ~ ~ 2.3C7V3 . c is in the range 0.08 to 0.2 giving 
a range of values for A from 7 to 19 for a real proteins in 
the laboratory, suggesting that droplet renormalization 
is not significant for the majority of cases. For the simu- 
lation models, on the other hand, the calculation gives A 
in the range 1.6 to 4 indicating that droplet renormaliza- 
tion is likely to be very significant for off-lattice models 
of random heteropolymers. 

TABLE I: Estimates of the reduced surface tension parameter 
A for some models studied by simulation. A 3p is the sharp in- 
terface approximation, while A% and A T are the results of ex- 
tracting A from the simulation data (using the RFOT theory) 
on the correlation length and relaxation times respectively. 

model k A 3p A^ A T 

p = 4 M = 3 @, J§ 1/6 10 4.8 - 
p = 4 plaquette Jo) 1/3 10 - 3.3 
q = 10 Potts jjj 1/6 8 

Computer simulations have been performed on a few 
spin-glass systems that exhibit a one-step RSB transi- 
tion in the mean-field limit including p = 4 spin models 
on different lattices j2^, ||(], and q — 10 state Potts 
glass [[l9| . Unlike the p = 4 spin models the Potts model 
did not display the characteristic signs of a glass transi- 
tion. An estimate for A for these systems based on the 
sharp short-range interface argument is clouded by the 
nearness of the mean-field Ta and Tk for these models. 
If we nevertheless assume a sharp planar interface, and 
take the free energy cost of forming the interface to be 
approximately one half of the energetic cost of random- 
izing the interactions across the interface, we would find 
rl<j = k\E k \ (or A = Airk \E K \ /T K ) with k the frac- 
tion of interactions per site that are broken in forming 
the interface and Ek the energy per site at Tk. The 
resulting estimates for A using £k and Tk from simula- 
tions (or from mean-field theory in the case of the Potts 
model) are shown in the table. Taken at face value, these 
suggest that the p = 4 spin models have a surface ten- 
sion only just sufficient to avoid substantial corrections to 
mean-field theory, while there may be significant renor- 
malization for the 10-state Potts model. Alternatively 
an estimate of the surface tension can be obtained using 
the observed growth in correlation length j2{| and relax- 
ation times [p0| which indeed have the functional forms 
predicted by RFOT theory. In doing this it is neces- 
sary to assume the theoretical results apply outside the 
region Tk < T < Ta. These estimates, also given in 
the table are much lower. These low values are consis- 
tent with the observation that Ta (where surface tension 
vanishes in finite ranged models) lies only slightly above 
Tk for these models in the mean-field limit; in particular 
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T A /T K = 1.009 for the 10-state Potts glass ||. We con- 
clude that the apparent lack of a glass transition in the 
10-state Potts model is due to a low effective microscale 
surface tension. Since |£"k| /?k = 25*o (in a random en- 
ergy model, with So the total entropy per site) it can be 
anticipated that A will increase approximately logarith- 
mically with the number of states in the Potts model but 
proportional to p in a p-spin model [ j32| . Achieving spin 
systems with a surface tension similar to that in a struc- 
tural glass would require a p w 20 (M = p—1) spin model, 
but a q — 10 3 (or even higher) state Potts glass. Unfor- 
tunately spin models that are relatively easily simulated 



appear to have substantial corrections to the current form 
of RFOT theory while real structural glasses and protein- 
like random heteropolymers have much smaller correc- 
tions. 
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